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ABSTRACT 

We report on the z-dropouts detected in two exceptionally deep ACS fields (-B435, V6O6) 2775 j and zs^o 
with 10a limits of 28.8, 29.0, 28.5, and 27.8, respectively) taken in parallel with the UDF NICMOS 
observations. Using an i — z > 1.4 cut, we find 30 i-dropouts over 21 arcmin 2 down to -zsso.ab = 28.1, 
or 1.4 i-dropouts arcmin -2 , with significant field-to-field variation (as expected from cosmic variance). 
This extends i-dropout searches some ^0.9 m further down the luminosity function than was possible 
in the GOODS field, netting a ^7x increase in surface density. An estimate of the size evolution for 
UV bright objects is obtained by comparing the composite radial flux profile of the bright i-dropouts 
(^850, ab < 27.2) with scaled versions of the HDF-N + HDF-S [/-dropouts. The best-fit is found with a 

(1+z) -0-53 scaling in size (for fixed luminosity), extending lower redshift (1 < z < 5) trends to z ~ 6. 
Adopting this scaling and the brighter i-dropouts from both GOODS fields, we make incompleteness 
estimates and construct az~6 LF in the rest-frame continuum UV (~ 1350A) over a 3.5 magnitude 
baseline, finding a shape consistent with that found at lower redshift. To evaluate the evolution in the 
LF from z ~ 3.8, we make comparisons against different scalings of a lower redshift B-dropout sample. 
Though a strong degeneracy is found between luminosity and density evolution, our best-fit model scales 
as (1 + z)~ 2 ' 8 in number and (1 + z) 0A in luminosity suggesting a rest-frame continuum UV luminosity 
density at z ~ 6 which is just 0. 38^007 x that at z ~ 3.8. Our inclusion of size evolution makes the 
present estimate lower than previous z ~ 6 estimates. 



Subject headings: galaxies: evolution 
function 



galaxies: high-redshift — galaxies: luminosity function, mass 



1. INTRODUCTION 

The unique deep z-band capabilities of the Hubble Space 
Telescope (HST) Advanced Camera for Surveys have un- 
covered a substantial population of z ~ 6 dropouts (e.g., 
Yan, Windhorst, & Cohen 2003). Early estimates indicate 
that the star formation rate at z ~ 6 is only somewhat 
lower than (only 0.17-0.77x that) at z ~ 3 (Bouwens et 
al. 2003b; Stanway, Bunker, & McMahon 2003; Dickin- 
son et al. 2004; Stanway et al. 2004; Bouwens et al. 2004, 
hereinafter, B04). Unfortunately, these early studies have 
all been on fields of somewhat intermediate depth, with 
strong incompleteness corrections beyond Zs5o,AB — 26.5. 
Since a simple extrapolation of the Steidel et al. (1999) 
z ~ 3 luminosity function to z ~ 6 predicts an of 
Z850,ab = 26, these samples are becoming incomplete just 
as the j-dropouts are starting to become prominent. More- 
over, the substantial biases against larger, lower surface 
brightness objects make it difficult to derive the size dis- 
tribution at z ~ 6 (e.g., B04 find that only ^25% of the 
brighter (7-dropouts make it into their z ~ 6 Great Ob- 
servatories Origins Deep Survey (GOODS) sample). It 
is therefore important to examine i-dropouts in a much 
deeper field. 



Here we report on observations of i-dropouts from two 
exceptionally deep ACS Wide Field Camera (WFC) paral- 
lels (hereinafter, UDF-Ps) that were taken in conjunction 
with Near Infrared Camera and Multi-Object Spectrome- 
ter (NICMOS) observations of the Ultra Deep Field (UDF) 
(GO-9803; Thompson et al. 2004). 9, 9, 18, and 27 orbits 
were taken in the F435W, F606W, F775W, and F850LP 
passbands (hereafter called -B435, Veoe, *775> an d Z85o) for 
each field over 9 pointings (the pointings were arranged in 
a 3 x 3 grid, each separated by 45"). These observations 
extend i-dropout searches some 0.9™ magnitude deeper 
than GOODS and to within 0.9 m magnitudes of the UDF, 
giving a preview of the results at these depths. We assume 
(n M ,Sl A ,h) = (0.3,0.7,0.7) (Bennett et al. 2003). 

2. OBSERVATIONS 

These observations were retrieved from STScI and sub- 
ject to the standard CALACS on-the-fly reprocessing (OTFR). 
Additional processing, including cosmic-ray rejection, align- 
ment, and final image combination, was performed by the 
"Apsis" data reduction software (Blakeslee et al. 2003). 
The latest photometric zero points were applied (Sirianni 
et al. 2004) along with a correction for galactic absorption 
(Schlegel, Finkbeiner, & Davis 1998). The 10cr limiting 
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magnitudes (for the deepest parts of the image) were ap- 
proximately 28.8, 29.0, 28.5, and 27.8 in the B A35 , V 60 6, 
Z775, and Z85o bands, respectively, with PSF FWHMs of 
0.08-0.09". 

3. ANALYSIS 

Areas whose V, i, and z exposures were shorter than 5, 
10, and 15 orbits were excluded from the analysis. Object 
detection and photometry was performed on the remaining 
area for each parallel (21 arcmin 2 , in total) using SExtrac- 
tor (Bertin & Arnouts 1996) with smaller scaled apertures 
(MAG_AUTO) to measure colors and larger scaled aper- 
tures to correct to total flux. For our i-dropout selection, 
we apply a similar i-dropout selection criteria to that used 
by B04 on the GOODS data set. We require objects be 
a 5.5 cr-detection within a 0.25"-diameter aperture, have 
a SExtractor stellarity less than 0.85 (non-stellar to high 
confidence), have (1775 — z S50 )ab colors redder than 1.4, 
and be a non-detection (< la) in the Vf306 band (within 
the smaller-scaled apertures). Several studies have demon- 
strated that a simple i — z cut can be an effective means of 
selecting z > 5.5 galaxies (Stanway et al. 2003; Bouwens 
et al. 2003b; B04; Dickinson et al. 2004). Since the surface 
densities of all known interlopers are expected to increase 
less quickly than bona-fide i-dropouts, we adopt the 11% 
found for brighter i-dropouts from GOODS (B04) as an 
upper limit to the contamination rate, with no spurious 
detections (running through a similar selection procedure 
on the negative images, no i-dropouts were found.) Given 
the similar incompleteness expected due to blending with 
foreground objects (verified through simulations, see B04), 
no net adjustment is made to the surface density. 

(a) Surface Density. With the above selection criteria, 
we recovered 30 i-dropouts to a limiting zsso.ab magni- 
tude of 28.1, 20 in the first parallel and 10 in the second, 
suggesting a significant degree of clustering, but not in- 
consistent with the cosmic variance (50% RMS) expected 
for fields of this size (Somerville et al. 2004; B04). The 
luminosities of i-dropouts in our sample range from 1.1 L* 
{z S 50.AB = 25.9) to 0.15 (z 8 50,ab = 28.1) (using the 
Steidel et al. (1999) value for at z ~ 3). Half-light radii 
(uncorrected for PSF effects) range from 0.1" to 0.3", or 
0.6 kpc to 1.7 kpc for the adopted set of cosmological pa- 
rameters. The positions, magnitudes, sizes, i — z colors, 
and SExtractor stellarity parameters of these objects are 
tabulated in Table 1. 

(b) Size Evolution. Cosmic surface brightness dimming 
results in significant incompleteness at z ~ 6 (Bouwens 
et al. 2003b; Giavalisco et al. 2004; B04). To treat this, 
it is necessary to have a good estimate of the size (sur- 
face brightness) distribution (independent of selection ef- 
fects). An empirical way of estimating this distribution 
is to model the observed sizes in terms of a scaled ver- 
sion of a lower redshift population and therefore derive 
the size evolution. For our low redshift comparison sam- 
ple, we adopt the z ~ 2.5 HDF-N + HDF-S [/-dropout 
sample (B04). The large Alog(l + z) increases our sensi- 
tivity to size changes. To artifically redshift this sample 
to z rsj 6-weighting by 1/V ma x and accounting for pixel k- 
corrections, surface brightness dimming, and detailed im- 
age properties like the noise and the PSF (as well as our 
selection procedure given above)-we employ the now well- 



established cloning machinery (Bouwens, Broadhurst, & 
Silk 1998a,b; Bouwens et al. 2003a,b; B04). Finally, this 
redshifted sample is compared with a sample of i-dropouts 
(^850, ab < 27.2) extracted from the ACS images degraded 
to the PSF of our [/-dropout images projected to z ~ 6 
(0.14" FWHM) with a S/N equivalent to 5 orbit V 60 6, 
10 orbit 2775, and 15-orbit z S50 observations (0.6 m deeper 
than GOODS: see sample #2 from Table 1 for the objects 
used here). 
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Fig. 1. — The mean radial flux profile for the 10 brightest 
z 850, AB < 27.2 i-dropouts (histogram, with ±lcr uncertainties) in 
our sample versus the predictions based upon different scalings of 
the z ~ 2.5 HDF-N + HDF-S [/-dropout sample (the 4775 images 
available for this sample). The shaded violet, green, and dark violet 
regions represent predictions based upon the (1 + z)°, (1 + z)~ 1,5 , 
and (1 + z)~ 3 size scalings, respectively. Clearly, the (1 + z)~ 1,5 
scaling is a much better fit to the observations than the other two 
scalings. This is quantified in the figure inset, which shows the 
expected la scatter [diagonal violet band) in the mean half-light 
radius vs. the (1 + z)~ m size scaling exponent m. Since the half- 
light radius for the mean profile is 0.85 ±0.06 kpc (plotted as a gray 

_ y 57 +o.53 

vertical band), the inferred scaling is (1 + z) ' -o.bq (shown as 
the solid black circle with la error bars) or 2.9^J |x from z ~ 2.5 
to z ~ 5.9 (the mean redshifts of the U and i-dropout samples, 
B04). This is in good agreement with the allowed range of scalings 
(m = 0.8 — 2.0) inferred from a recent analysis of the GOODS data 
(B04) (shaded in light orange). 

We consider size scalings of the form (1 + z)~ m (m = 
— 3) (for fixed luminosity) and quantify the differences 
in terms of the mean radial flux profile, which provides 
a measure of the light within different circular annuli (see 
B04 for the computational details). The results of the pro- 
cedure are shown in Figure 1 for the (1 + z)°, (1 + z) -1 - 5 , 
and (1 + z) -3 scalings, and it is obvious that while the 
(1 + z) -1 - 5 scaling provides a good fit, the (1 + z)° and 
(1 + z) -3 scalings produce profiles which are too broad 
and narrow, respectively (see Figure 2 for a visual compar- 
ison of the brighter dropouts with the no-evolution m = 
scalings). To use these results to place confidence limits 
on the allowed scalings, we computed the mean size and 
variance for a 10-object sample of intermediate magnitude 
i-dropouts. The result is shown in the inset to Figure 1 
as the shaded violet region. Given the observed half-light 
radius 0.85 ± 0.06 kpc (corrected for PSF effects) (shaded 
grey region), we infer an m of 1.57^q'5q (shown in the in- 
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Fig. 2. — Postage stamp images (zsso-band) of 10 randomly-selected i-dropouts from the UDF-Ps (zsso.AB < 27.2) versus similarly selected 
no-evolution projections of our lower rcdshift HDF [/-dropout sample (B04). The observed i-dropouts are both smaller and of higher surface 
brightness than the cloned [/-dropout sample (a 2.9^q'|x evolution in size, see Figure 1). 



set as the solid circle with error bars), similar to what was 
found from our GOODS study (light orange region) where 
m was found to be between 0.8 and 2.0 (B04). 

(c) Luminosity /Density Evolution. To interpret the i- 
dropout counts and therefore estimate the evolution of the 
LF, we use different projections of the lower redshift B- 
dropout sample from the GOODS fields (B04). We explore 
scalings in both luminosity (e.g., (1+z) 1 ) and number (e.g., 
(1 + z) n ), assessing these scalings by comparing against 
both the UDF-Ps (degraded to the B-dropout PSF pro- 
jected to z ~ 6, 0.1" FWHM with a S/N equivalent to 5 
orbit V(306, 10 orbit 1775, and 15 orbit zsso observations: 
see sample #3 from Table 1 for the objects used here) 
and the GOODS fields (B04) (corrected up to the same 
completeness level as that of the UDF-Ps), which provide 
important constraints at bright magnitudes. As per the 
results of the previous subsection, sizes of individual B- 
dropouts are scaled as (1 + z) -15 (increasing their surface 
brightness). Interpolating between the simulation results 
and calculating a reduced x 2 for each comparison, the 68% 
and 95% likelihood contours can be determined, taking as 
inputs the GOODS z-dropouts alone (shown in the Fig- 
ure 3 inset as the green contours) and GOODS combined 
with the UDF-Ps (black contours). While the results from 
the UDF-Ps are helpful in breaking the degeneracy be- 
tween luminosity and density evolution, small differences 
in normalization between the bright (from GOODS) and 
faint (from the UDF-Ps) counts could affect the results 
(due to cosmic variance). Deep wide- field images will 
be needed to ultimately establish the form of the evolu- 
tion. The (1 + z)~ 2 - 8 scaling in number and (1 + z)° 
scaling in luminosity preferred in this analysis point to- 
wards az~6 luminosity density (~ 1350A) that is just 
0.38+n 07 x that at z ~ 3.8 (error bars computed from 2D 
likelihood contours). Applying this factor to the z ~ 3.8 
rest-frame continuum UV luminosity density determined 
by B04 and using standard assumptions to convert this 
into a star formation rate density (e.g., Madau, Pozzetti, 
& Dickinson 1998), the star formation history integrated 
down to 0.2 L* can be plotted and compared with several 
previous determinations (Figure 4). The present result is 
noticeably lower. This is almost entirely the effect of size 
(surface brightness) evolution (see Figure 1) on the esti- 
mated completeness of the z ~ 6 census (and not due to 
our greater depth). Without size evolution from z ~ 3.8, 
the completeness correction and inferred luminosity den- 
sity are some ~ 1.9 x higher than the current result (closer 
to previously published estimates: Bouwens et al. 2003b; 
Giavalisco et al. 2004; B04-see Figure 4). 
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Fig. 3. — The observed number of i-dropouts per half magnitude 
interval observed in the UDF-Ps (solid black histogram shows the 
counts from the data degraded to match the projected PSF of the 
lower redshift samples [the undegraded data have an additional 8 
objects beyond 27.5, see Table 1]) versus that predicted by a lumi- 
nosity (blue-shaded region) and density (red-shaded region) scaling 
of the lower redshift B-dropout population. Each was scaled in size 
by (1 + z)~ 1,5 to match the evolution found in Figure 1 (the size 
scaling increased the B-dropout count predictions by 30-100%). The 
green histogram shows the number of i-dropouts per half magnitude 
interval observed in the GOODS fields using similar selection cri- 
teria (B04) corrected up to the completeness level of the UDF-Ps. 
The figure inset shows the 68% and 95% confidence intervals for lu- 
minosity evolution (parametrized as (l + z) 1 ) and density evolution 
(parametrized as (1 + z) n ) taking as inputs the GOODS data alone 
(green contours) and the GOODS plus UDF parallel data (black 
contours). While the observations (bright i-dropouts from GOODS 
and faint i-dropouts from the UDF-Ps) prefer a model where I = 0.1 
and n = —2.8, pure luminosity evolution or pure density evolution 
cannot be excluded (at 95% confidence). 



(d) Luminosity Function. It is also useful to construct a 
simple estimate of the LF to examine the shape and look 
for evolution in the faint end slope. For our procedure, we 
follow Bouwens et al. (2003b) and treat the determination 
of the luminosity function 4>(M) as a deconvolution prob- 
lem on the integral J (f>(M(m, z))p(m 7 z)^dz — N(m) 
where pirn, z) is the selection efficiency (including vari- 
ous factors ranging from galaxy color to photometric scat- 
ter to detectability) , N(m) is the number counts, and the 
absolute magnitude M here is given as a function of the 
apparent magnitude m and the redshift z (using an av- 
erage color dropout to perform the fc-correction). Adopt- 
ing N(m) from the primary data sets and taking p(m, z) 
from the i?-dropout cloning simulations (scaled in size as 
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Fig. 4. — A history of the star formation rate density assuming 
no extinction correction, integrated down to 0.2L*. The red point 
at z ~ 5.9 represents our estimate based upon the i-dropouts found 
in both GOODS and the UDF-Ps (§3(c)). A Salpeter IMF is used 
to convert the luminosity density into a star formation rate (e.g., 
Madau et al. 1998). Comparison is made with the previous high 
redshift determinations of Lilly et al. (1996) (open squares), Steidel 
et al. (1999) (green crosses), Giavalisco et al. (2004) (solid black 
diamonds), and B04 (solid magenta circles). The top horizontal axis 
provides the corresponding age of the universe. The uncertainties for 
all determinations are formal errors. The inclusion of size evolution 
makes the present estimates lower than previous z ~ 6 estimates. 



(1 + z) -15 ) (uncertain on the 20% level due to uncertain- 
ties in the size distribution), we are able to compute the 
LF and provide a simple fit to the Schechter function (Fig- 
ure 5). Note that the relative normalization of the bright 
end of the LF (Mi 350 .ab < —19-7), determined primarily 
from GOODS, and the faint end of the LF, derived from 
the UDF-Ps, could vary by 36% RMS (the expected cos- 
mic variance for a CDM-type power spectrum normalized 
to the high redshift observations, see B04), affecting the 
shape of the LF. 
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Fig. 5.— The rest-frame continuum UV (~ 1350A) LF at z ~ 6 
determined from the GOODS fields (at M1350, ab < —19-7) and the 
UDF-Ps (red circles with ler error bars), with the best-fit LF to a 
Schechter function also shown (M1350 ab,* = —20.26, rf>* = 0.00173 
Mpc -3 , and a = -1.15). The Steidel et al. (1999) z ~ 3 LF (dot- 
ted green line) is presented for comparison. A fit with a = —1.6 (as 
observed at z ~ 3, Steidel et al. 1999) is also consistent (at 68% con- 
fidence) with the data, and yields values for Mi35o,AB,» (=—20.77) 
and <j>* (=0.00098 Mpc -3 ) more suggestive of the interpretation 
given in §3(c) (an evolution in density). 



4. DISCUSSION 

The additional 0.9 m we have in zsbo depth (and l.l m 
in i 775 ) over the GOODS field allows us to probe further 
down the LF and obtain a wider distribution of sizes at 
bright magnitudes than was possible from previous data. 
This depth largely compensates for the (1 + z) 4 cosmic sur- 
face brightness dimming, allowing us to set a more direct 
constraint on object sizes at z ~ 6 (cf. B04). The results 
here are consistent with the (l + z) -1 ' 5 scaling of size found 
at lower redshift (1 < z < 5) for objects of fixed luminos- 
ity (Ferguson et al. 2004), extending the trend to z ~ 6. 
This evolution in size has important implications for the 
luminosity density inferred at this epoch, suggesting that 
the star formation rate density has increased significantly 
from z ~ 6. 

In this work, we present the first robust luminosity func- 
tion at z ~ 6 and find no radical difference relative to lower 
redshift. While the UDF will add depth and further con- 
straints, the present fields provide a crucial assessment of 
the cosmic variance. 
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Jon McCann for their assistance and the referee for valu- 
able comments. We acknowledge the support of NASA 
grant NAG5-7697 and NASA grant HST-GO09803.05-A. 
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Table 1 

z ~ 6 objects from the udf-ps. a 



Object ID Right Ascension Declination Z850 i — z S/G r hl(") Sample 



UDFP1-4109-2788 


03: 


:32: 


:42 


.692 


-27: 


56:55 


.35 


25 


,9±0, 


.1 


1 


.8 





.00 


0. 


.34 


3 


UDFP1-3629-862 


03: 


:32: 


:49 


.638 


-27: 


56:27 


.47 


26 


,2±0, 


.2 


2 


.0 





.02 





.26 


3 


UDFP1-3851-2438 


03: 


:32: 


:43 


.958 


-27: 


56:43 


.86 


26 


,3±0, 


.1 


1 


.6 





.18 





.11 


1,2,3 


UDFP1-4650-3354 


03: 


:32: 


:40 


.706 


-27: 


57:24 


,28 


26 


,4±0 


1 


1 


.7 





.02 


0. 


18 


1,2,3 


UDFP1-1546-1408 


03 


:32: 


:47 


047 


-27: 


54:48 


.09 


26 


,9±0 


.1 


1 


.4 





.00 





19 


2 


UDFP1-2954-1152 


03' 


:32: 


:48 


.367 


-27: 


55:54 


.81 


26 


,9±0, 


.2 


1 


.5 





.37 


0. 


.11 


2 


UDFP1-4040-2514 


03: 


:32: 


:43 


.709 


-27: 


56:50 


.93 


27 


,0±0, 


.2 


1 


.9 





.01 


0. 


19 


1,2,3 


UDFP1-2309-1628 


03: 


:32: 


:46 


.409 


-27: 


55:24 


.31 


27 


,1±0, 


.1 


>2 


.5 





.02 


0. 


.11 


1,2,3 


UDFP1-4784-3382 


03: 


:32 


:40 


.634 


-27: 


57:31 


.07 


27 


,2±0, 


,2 


1 


.6 





.02 


0. 


.14 


1,2,3 


UDFP1-3350-3425 


03 


:32: 


:40 


111 


-27: 


56:22 


.32 


27 


,4±0. 


.2 


1 


.5 





.01 





12 


1,2,3 


UDFP1-4413-3612 


03: 


:32 


:39 


671 


-27: 


57:13 


.38 


27 


,4±0. 


.2 


>2 


.1 





.02 


0. 


.11 


1,2,3 


UDFP1-3144-2410 


03 


:32: 


:43 


.876 


-27: 


56:08 


.48 


27 


,4±0. 


.2 


>1 


.8 





.01 





19 
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a Selection criteria: (1775 — ^85o)as > 1-4 plus null detection (< 2a) in V606- All limits are 2<r. "S/G" denotes the SExtractor 
stellarity parameter, for which = extended object and 1 = point source. Several different samples were used in the analysis. 
Sample #1 was generated by applying our selection criteria to the undegraded images. Sample #2 and #3 were generated by 
applying these selection criteria to the images degraded to match the PSF of the U and B-dropout samples, respectively, projected 
to z ~ 6 (0.14" FWHM and 0.10" FWHM). Images from the latter samples were degraded to a uniform S/N equivalent to 5 
orbit Vaofi, 10 orbit 4775, and 15 orbit 2850 observations. The union of all three samples (38 objects) is larger than the undegraded 
sample (sample #1) due to object blending and photometric scatter. UDFP1 and UDFP2 refer to the first and second UDF 
parallels, respectively. 



